Bistable Spin Currents from Quantum Dots Embedded in a Microcavity 
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We examine the spin current generated by quantum dots embedded in an optical microcavity. 
The dots are connected to leads, which allow electrons to tunnel into and out of the dot. The spin 
current is generated by spin flip transitions induced by a quantized electromagnetic field inside the 
cavity with one of the Zeeman states lying below the Fermi level of the leads and the other above. 
In the limit of strong Coulomb blockade, this model is analogous to the Jaynes-Cummings model 
in quantum optics. We find that the cavity field amplitude and the spin current exhibit bistability 
as a function of the laser amplitude, which is driving the cavity mode. Even in the limit of a single 
dot, the spin current and the Q-distribution of the cavity field have a bimodal structure. 

PACS numbers: 42.50.Pq,73.63.Kv,78.67.Hc 
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I. INTRODUCTION 

Optical bistability (OB) occurs when a nonlinear di- 
electric is placed inside of optical resonator, which pro- 
vides a feedback mechanism for the light. The output 
intensity from the resonator can, as a result, have two or 
more stable values for a given input intensity that is driv- 
ing the cavity and can be made to switch between these 
two outputs by varying the input intensity beyond the 
bistable regioni. The interest in OB systems started in 
the 1960's with Szoke's et al theory of absorptive optical 
bistability where the bistability is a result of absorption 
by the dielectric medium 2 -. OB was first observed ex- 
perimentally and explained by Gibbs et al? for a cavity 
containing a medium with a nonlinear index of refrac- 
tion and no absorption or gain (dispersive bistability). 
Interest in OB has been stimulated by its practical ap- 
plications including optical switches, logic gates, and op- 
tically bistable memory devices^* 6 -! 7 -! 8 -. In addition to 
device applications, OB is also interesting because it ex- 
hibits novel physical behavior such as phase transitions 
between stationary but non-equilibrium states^ 1 ^. 

In a completely independent development, spintron- 
ics has emerged as new field in which the spin degrees 
of freedom of charge carriers in solid state devices are 
exploited for the purpose of information processing. Ma- 
nipulation of the spin degrees of freedom rather than 
the charge has the advantage of longer coherence and 
relaxation times 11 . In order to manipulate the spin de- 
gree of freedom for the purpose of information process- 
ing, there is a demand for efficient and readily fabricated 
spin devices such as spin batteries, spin filters, spin tran- 
sistors, etc.. Much of this work has focused on ways 
to generate pure spin currents in semiconductor nanos- 
tructures using, for example, the extrinsic Spin-Orbit 
(SO) interaction^ 3 -, Rashba SO interactions^, optical 
absorption^, Raman scattering^, shape deformations of 
open quantum dot a 17 ' 18 , as well as various types of quan- 



tum pumps 
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Pure spin currents, I s = s(Jf — Jj_), 



are the result of an equal number of spin up (f ) and spin 
down {[) charge carriers moving in the opposite direction 



so that the charge current is zero, I c = q{I^ + IjJ. Here, 
I a are the spin polarized particle currents, s = H/2 the 
spin of the particle, and q = e the charge. 

Another model for a spin battery that has recently 
been proposed is electron spin resonance (ESR) in a 
quantum dot connected to leads, which generates a pure 
spin current by spin flip transitions when there is a large 
Zeeman splittin g 23 ' 24 . A classical transverse magnetic 
field was used to induce the spin flips and hence cre- 
ate the spin current. We recently extended this model 
by considering spin flips induced by a quantized mode 
of an optical microcavity 2 ^,. In this two-photon 
Raman transition via an intermediate charged exciton 
(trion) state was used to induce spin flips. The spin cur- 
rent was found to be significantly larger in our case than 
for a classical undepleted field as a result of the cavity 
decay. Also, the shot noise exhibited a rich structure that 
was consequence of the discrete photon numbers in the 
cavity. 

In our previous work, a spin flip from the lower to 
upper Zeeman states involved the absorption of a photon 
from a classical pump laser and creation of photon in the 
cavity mode. As a result, the cavity field was built up out 
of the vacuum without any need for driving the cavity^. 
Here, we consider the reverse process in which a photon 
from the cavity mode must be absorbed in order to flip 
the electron spin. This requires that the cavity mode now 
be driven by an external source. 

We show here that the driven cavity system exhibits 
absorptive OB for the amplitude of the cavity field. Be- 
cause the spin current is a function of the cavity field 
amplitude, the spin current also exhibits bistability as 
function of the amplitude of the driving field. We study 
the limit of both a single quantum dot coupled to the 
cavity as well as N 3> 1 quantum dots coupled to the 
cavity. In the later case, we can use a semiclassical treat- 
ment for the dots and cavity field 26 . This behavior in- 
dicates that the system could be used as an optically 
controlled spin current switch. While optical bistability 
has been studied in quantum wells embedded in semicon- 
ductor microcavities^ 7 -, this is the first study of OB in the 
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presence of electrical transport through quantum dots. 

In Section II, wc outline our model. In Section III, we 
analyze the steady state solution of the master equation 
for a single quantum dot coupled to the cavity. In Section 
IV, we consider the semiclassical solution for N 3> 1 dots 
coupled to the cavity mode. We examine the effect that 
inhomogeneous broadening of the dots' Zeeman splitting 
and variations of dots' vacuum Rabi frequency have on 
the bistability. 

II. MODEL 

We consider self-assembled quantum dots embedded 
in a high-Q microcavity, as depicted in Fig. 1. Strong 
coupling between individual self-assembled and interface 
fluctuation quantum dots with a single mode of an op- 
tical microcavity has recently been achieve d 28 ' 29 . Here 
we are interested in simultaneous coupling of dots to a 
cavity mode and electrical transport through the dots 
due to tunnelling from a doped reservoir. Although 
self- assembled quantum dots are usually used for opti- 
cal studies, there have been several experimental studies 
of electrical transport through individual self-assembled 
InAs quantum dot a 30 ' 31 ' 32 as well as through a thin sheet 
containing a large number of InAs dots 33 . Along a simi- 
lar line, the ability to control the tunnelling of electrons 
or holes between self-assembled dots and a doped GaAs 
reservoir by a gate voltage combined with simultaneous 
spectroscopic studies of these charged quantum dots has 
been demonstrate d 34 ' 35 . 

We assume that there are two electron reservoirs at 
chemical potential, /i^ = = /1, that are coupled to the 
dots via tunnelling. We note, however, that our results 
are equally valid in the case of only a single reservoir. 
Only a single empty orbital energy level, £j, of the dot 
lies close to /1. Here, the subscript i denotes the particu- 
lar dot. The Zeeman splitting between the two electron 
spin states is Aj = e^j — e^j = g x (J^bB where B is a 
static magnetic field along the x-axis that is perpendicu- 
lar to the growth direction (z). [Ib is the Bohr magneton 

and g$ is the electronic g-factor of the i th dot along the 
direction of the magnetic field. The energy levels satisfy 
e.j.j = £j — Aj/2 < /i < e-jj = Ei + Aj/2 so that only 
spin up electrons can tunnel into the dot and only spin 
down electrons can tunnel out. In the limit of very large 
Coulomb blockade energy, only a single electron from the 
reservoir can occupy the dot resulting in the bare Hamil- 
tonian for the dots and cavity field, 

N 

i=i 

where c iiCT (cJ CT ) are annihilation (creation) operators for 
electrons in dot i with spin a in the x-direction of the 
magnetic field. 

Transitions between the different electronic Zeeman 
states of a dot are induced via a two-photon Raman tran- 
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FIG. 1: (a) Raman transition between the dot Zeeman states, 
I t, I), via an intermediate trion state, |+t), induced by a laser 
with frequency lul and a cavity mode with frequency u> C av- 
The spin eigenstates along the direction of the magnetic field 
are superpositions of spin eigenstates in the growth direction, 
c' ±z |0). (b) Schematic of a single quantum dot indicating Zee- 
man energy levels in the dot and allowed tunnelling between 
leads and dot. 

sition involving a strong laser field that may be treated 
classically and a quantized mode of the microcavit y 25 ' 36 . 
The two optical fields couple the electron spin states to 
a higher energy charged exciton state (trion) 35 ' 37 ' 38 ' 39 . 
The lowest energy trion states excited by cr + and a~ cir- 
cularly polarized light are | + U) = cj ^.cj ^ hi +3 / 2 |0) an d 
I — ti) = cj jet Ji\ _ 3 / 2 |0) , respectively. They consist of 

an electron singlet state and a heavy hole, where h\ ±3 ^ 2 
are heavy hole creation operators with spin projections 
±Ti3/2 along the z-axis and |0) is the empty dot state. 
The a + polarized pump laser with frequency u>i and Rabi 
frequency fi; couples each of the electron spin states to 
the I +ti) trion state. Similarly, we assume that the cav- 
ity field, with vacuum Rabi frequency g ca v,i for each dot 
and frequency uj cav , is also <j + polarized due to either 
the cavity construction or because it is driven by a er + 
pump as discussed below 4 ^. When the two fields are far 
detuned from the creation energy for the | + ti) state, the 
intermediate trion state can be adiabatically eliminated 
to give 

N 

Hj = inJ29i(ac[iCi^ iuit ~ h.c). (2) 

i=l 

where — g C av,i^l/^^-R an d A/j is the detuning from 
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the trion state. In deriving Eq. [2l we have assumed that 
W - {u C av < W + {u ca v ~ ^l)\ so that the non- 

resonant terms d^ct , di ^e~ tuJlt + h.c. can be neglected. 
Since electrons enter the dot in the spin | state, a pho- 
ton must be absorbed from the cavity mode in order to 
generate a spin current. It is therefore necessary to pump 
the cavity field. We assume that the cavity is driven by 
a classical source oscillating at frequency uj p , 



ihe(a't 



h.c), 



(3) 



which generates a coherent state in an empty 
cavity2^£. 



The explicit time dependence can be removed from 
Hj and Hp by transforming to a rotating frame for 
the field operators, a ~ Ae~ luJpt and the dot operators, 
Cj,T = Ci,| exp(-i(ui - u p )t/2) and Cjj = C hi exp(i(uji - 
uip)t/2). The Hamiltonian in the rotating frame is H' = 
Hq + Hp + Hj , 



N 



K = fi(u cav - lo p )A^A + (ei(C{ T C ilT + Cl L C itl ) + (A, + w , - w p )(Cj, - Cj T C ilT )/2) (4) 

i=l 

N 

H'j + H'p = ih^2g i (Acl l C i ,- [ -h.c.) + ihe(A^ -h.c.) (5) 



r 



From Eq. 0] one can clearly see that the resonance con- 



ditions are uj c 



and A, 



u>i . We assume that 



former condition is always satisfied, while the latter con- 
dition cannot be satisfied for all dots due to variations in 
the magnetic moments between dots. 

The dynamics of the system can be described in terms 
of the density operator, p for the cavity plus dots. The 
master equation for p is given by, 



-i[H',p] 



,(A^Ap - 2Ap& + piU)/2 + p\i 



The first term describes coherent dynamics of the coupled 
QD-cavity system, the second term stands for the cavity 
decay^ 2 -, and the third term describes QD-lead coupling. 
Here we assume that the Coulomb blockade is so large 
that a second electron cannot tunnel into the dot if there 
is already one electron in the dot. The lead-dot coupling 
is most easily expressed in terms of the matrix elements 



of the density operator, p 



(n,m) 



where 



\<Ji,n) represents a state with n photons in the cavity 
and <ii = 0, |, I corresponding to no electrons, one spin 
up, or one spin down electron, respectively, on the i th 
dot. The specific form of the master equations for the 
lead coupling ar o 24 ' 25 



. (n,m) I 
P0i,0i llead 

. (n,m) I 
P'\ i \ i \lead 

. (n,m) I 
P I i I i \lead 

. (n,m) I 
P-K K \lead 



p(-)„( n .™) p(+) _(","») 
i "ii,U i POifii 

p(+) „(n,m) 
L i POtfii 

l i p u,u 
_r<~) n ( n ' m ) in 



(7) 
(8) 
(9) 
(10) 



Here, v\ ' = ^ + rj ^ is the rate at which spin down 
electrons tunnel out of the dots into the left and right 



^i + R ' s ^ ne ra ^ e a ^ w hich spin 
up electrons tunnel into the dots. (The subscripts L and 



leads and t[ +) = 



R denote the tunnelling rates for the left and right leads, 
respectively.) We also assume that the coupling between 
the left and right leads and the dots are the same and that 
the tunnelling between the leads and the dot is spin in- 

]?W = T. It is worth 



= r 



(-) 



,(+) 



dependent, x i h - x l L - x i R - x . fl 
pointing out that because Eqs. [7lfT0l take into account 
large Coulomb blockade in the dots, they are slightly dif- 
ferent from other master equations for dots that assume 
noninteracting electrons 4 ^. 



III. SINGLE QUANTUM DOT 

First we consider the limit of a single quantum dot 
and exact two-photon resonance, A = lu cciv — u>i. We nu- 
merically solve for the steady state density matrix by first 
expressing the density matrix in vector form, p^™'™" 1 — > p, 
and rewriting Eq. [S]in matrix form, 



dp I dt = M p 



(11) 



The steady state solution, p^ ^ , is given by the eigenvector 
of M with zero eigenvalue 4 ^. 

The steady state behavior of the system can be charac- 
terized by the Q-distribution for the intracavity optical 
field 1 -! 4 ^, Q{a) = J2a=o i i( a i a \P^\ a > a ) / n where |a) is 
a coherent state, a\a) = a\a). The advantage of the Q- 
distribution is that it is positive semi-definite and can be 
used to make comparisons to classical phase space prob- 
ability distributions. Figures 2 and 3 show Q(a) as well 
as the probability distribution for the photon number, 
P(n). One can clearly see two peaks corresponding to 
the two stationary average cavity field amplitudes. 

The average spin current, I s — s(7j — ij), is the same in 
the left and right leads with the stationary currents given 
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FIG. 2: Q-distribution vs. Re[a] and Im[a] for V ca v = 0.2r, 
g = 1.4T and e = 0.5r. 




FIG. 3: Photon probability distribution, P(n), vs. cavity field 
photon number, n, for same parameters as Fig. 2, T ca „ = 
0.2r, gr = 1.4T and e = 0.5r 
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are the 



populations of the dot after tracing over the state of the 
cavity held. The average spin current can be expressed 
in terms of expectation values of the cavity held using 
Eq. m 



(I s ) = 2s(2eRe[ 



(12) 



Since the bimodality in Q(a) is in the amplitude, \a\, 
both i?e[a] and ala have two most probable values corre- 
sponding to the two peak locations, cxj for j = 1, 2. Con- 
sequently, individual measurements of the spin current 
would give results clustered around two most probable 
values 



I S j = 2s(2eRe[aj] - T cav a*a 3 ). 



(13) 



By contrast, the ensemble averaged spin current of Eq. 
[His (I s ) = Tr[I sP {Q) } Pa P^i + ^2/5,2 where P j are 
the total integrated probabilities for the two peaks in 
Q(a). Because of the large quantum fluctuations associ- 
ated with a single dot, aw are not truly stable points of 



the system. Instead, this system should behave in a sim- 
ilar manner to individual atoms in optical cavities, which 
exhibit stochastic jumps between the two stationary cav- 
ity field values induced by quantum nois o 26 i 41 ' 45 > 46 . 



IV. SEMICLASICAL BISTABILITY 

Here we consider the equations of motion for the field 
and dot expectation values, 



a = (A) 

Si = (Cj^C'i^) 

P,i = (ct&i) 

Pi,] = (^.At), 



(14) 
(15) 
(16) 
(17) 



which can be derived from Eq. [5] using (O) — Tr[Op\. 
In order to have a hnite system of equations, we must 
factorize the expectation values involving the cavity field 

and dots, (i^Ai) ~> (^ f ) AAi> = a * s " Thc 
resulting equations are then, 



Pi, I 
Pi,1 

a 



gi(sia* + s*a) - Tp itl (18) 

-g l {s l a* + s*a) +T(l-p ld -pij) (19) 
JV 

-^giSi-T cav a/2 + e (20) 



i=i 



^ = -iSiSi + gia(pi^ - p t ^) - Ts l /2 (21) 

where 5i = Aj + lui — lu p . This 'semiclassical' factoriza- 
tion ansatz amounts to neglecting quantum mechanical 
correlations between the cavity held and dots and is usu- 
ally assumed to be valid in the limit of large 'classical' 
systems 2 -. Experiments with atoms in optical cavities in 
the strong coupling regime showed good agreement with 
such a semiclassical theory for TV > 15 atoms except for 
very close to the end points of the bistable region 47 . We 
therefore restrict our treatment to N ^> 1 dots. 

In order to derive analytic expressions, we will first 
neglect inhomogeneous broadening, Si = 0, and varia- 
tions of the vacuum Rabi frequency due to random po- 
sitions of the dots relative to the cavity mode, gt = go. 
We will account for these effects later by numerical av- 
eraging over Si and gi. We then define new variables 
for the total population and polarization, P a — J2iPi,<? 
and S = Si. By introducing the polar representation, 
S = \S\e i0 , a = \a\e lv and e = |e|e^, and using the 
positive dchnitcncss of P ai one sees that the phases are 
locked, <p = 4> = 9. 

The differential equations for \a\, \S\, and P„ have 
steady state solutions given by 



\2glY cav \a\] 



\S\i 



24 5o 2 |6|H 2 + (4.g 2 rW- 



2r 2 |e| 
2|e|- 



r 2 r 



a 



2ffo 



■j)\<*\j 
(22) 

(23) 
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Pa = 



-r ca >|?+2|el|a| 



P Tj = N-2P\ 



(24) 
(25) 

land the 



Here, \a\j for j — 1, 2, 3, are the roots of Eq. 
steady spin current is given by I S j = 2sTPij 

In the limit of negligibly small cavity decay, T cav <C 
|e|, T, g , Eq. (j2"2"]) becomes quadratic with the two 
stationary solutions, 



a 1,2 



tn 

12H 



1± 



N 2 g 2 



(26) 



0. 



for < \2\e\ 2 /N 2 gl < 1. By contrast, in the limit V 
the only nonzero solution for the cavity field is |a| = 
2|e|/r cal) . This limit is different from the case of atomic 
OB where bistability in the phase of a exists in the limit 
that the spontaneous emission rate goes to zero 45 ' 46 . The 
difference arises from the fact that T represents both the 
pumping rate and the decay rate for electrons in the dots. 
When r = 0, electrons cannot tunnel into the dots and 
interact with the cavity, which effectively decouples the 
dots from the cavity field. 
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FIG. 4: (a,c) Cavity field amplitude log 10 (|a|/JV) vs. \e\/N.(b, 
d) Spin current, I s /2sFN vs. \e\/N. For (a) and (b), g = 
1AT and T cav = 0.2r. For (c) and (d), g = 2.0F, Tcav = 
0.5F. In all cases the number of dots is N = 50 and \e\/N is 
plotted in units of V. 

Solutions of Eq. [2U are presented in Fig. 4 as a func- 
tion of the driving field, e. Three positive real solutions 
of Eq. [22] exist for e in the interval |e|i < |e| < | e 1 2 , where 



|e|i,2 = 




1 fN , 



No 
N 



3/2 



(27) 



where N = 2IT c(M ,/<?q i s the critical dot number. By 
analogy to atomic cavity QED, it is the number of dots 
necessary to significantly modify the resonant properties 
of the cavity 26 . The requirement that (1 — N /N) 3 / 2 



be real implies that Nq < N, which reduces to the re- 
quirement of strong coupling to the cavity when TV = 1, 
as one would expect. Eq. [27] corresponds to the range 
of I e I for which bistability occurs. One can clearly see 
from Fig. 4 that the bistability of a also leads to bista- 
bility and hysteresis for the spin current. We have an- 
alyzed the stability of the three steady state solutions 
as outlined in the appendix and determined that the 
two positive slope solutions in Fig. 4 (represented by 
the solid lines) are stable while the negative slope solu- 
tions (dotted lines) are unstable. For the cavity field, 
the lower stable branch corresponds to a nearly empty 
cavity, a ~ but with a finite increasing spin current, 
Is « up to I s m 2siVr x 0.2. Since g ,T > T cav , this 
solution corresponds to a photon entering the cavity and 
being almost instantaneously absorbed by a dot, which 
then leads to the creation of one unit of spin current. The 
upper stable branch corresponds to a large cavity field, 
a « O.liV up to a 10./V, with a constant spin current 
of Is — 2sTN/3. This value of the spin current cor- 
responds to a saturated transition in the dot with equal 
probabilities for the two spin states along with the empty 
dot state, p °^ 0i = P^i t = Puu = 1/3- I n this case, a 
spin current per dot of 2sr/3 was previously found to be 
the maximum spin current that could be generated by a 
single dot using a classical field 2 ^. 

In order to account for variations in Si and gi as a result 
of the random sizes and positions of the dots during the 
growth process, we assume that the variations can be 
modelled using Gaussian probability distributions, 



P d {5 l ) 
P g (9i) 



1 



<TdV 2tt 

1 



Og\/ 2tt 



(28) 
(29) 



The contribution of each of the dots to the cavity field 
appears as a summation over the dot polarizations in 
Eq. 1201 which must now be replaced with an inte- 
gration over the probability distributions, J^i 9i s i ~ * 
N J J dSidgiPd(Si)P g (gi)giSi. The resulting equation for 
the steady state cavity field is then 



|e|-r ca „|a|/2 = N I d8 t dg l Pd{5 l )P g {g l )f{5 11 g 1 \ , i0) 

s?l«|r 



f(Si,9i) 



6 ff 2 |a| 2 +r 2 /2 + 2Sf 



(31) 



The cavity field amplitude can then be used to calculate 
the spin current using I s = 2sTJ2iPi.l — 2s(2|a||e| — 
T cav \a\ 2 ). Note that in deriving Eq. [3D] 

we have made 

use of the fact that the sine of the phase difference be- 
tween the cavity and driving fields, sin(0 — (p), vanishes 
when averaged over a probability distribution that is even 
about Si = 0. One can see from Eq. [301I3T1 that for small 
cavity fields gi\a\ <§; T, the homogenous broadening of 
the dot levels due to the leads dominates over the inho- 
mogeneous broadening of the dot levels when <C T. 
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On the other hand for large cavity fields, gi\a\ ^> V, ad, 
power broadening dominates over homogeneous and in- 
homogeneous broadening of the dot levels and the cavity 
field is independent of ad- 
Figures 5 and 6 show the cavity field amplitude and 
spin current as a function of the driving field amplitude, 
e, for different ad and a g . For the Rabi frequencies, a g 
was chosen to be 20% or 40% of the mean while for the 
Zeeman splitting we chose ad > T for all curves. One can 
see that for increasing ad and a g the range of e where 
bistability occurs is reduced. This can be understood as 
the result of fewer dots that are both resonant with the 
cavity mode and strongly coupled to it. 
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FIG. 5: Cavity field amplitude, \ot\/N, vs. |e|/JV in units 
of T for TV = 50 dots, g = 1.6r, and F cav = 0.4F The 
different curves correspond to: ad = 2T and a g = 0.2go (solid 
line); ad — 2T and a g — OAgo (dashed line); ad = 4r and 
a g = 0.2<;o (dotted line); ad = 7T and a g = 0.2go (dashed dot 
line). 



switch between the two stable states by varying the 
driving field amplitude, e, beyond the endpoints of the 
bistable region^. One might then envision using this sys- 
tem as an optically controlled spin current switch, which 
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FIG. 6: Spin current, I s /2sTN, vs. |e|/JV in units of T for 
the same parameters as Fig. 5: TV = 50 dots and go = 1.6T, 
F C av ~ OAT. The different curves correspond to: ad = 2T 
and a g = 0.2go (solid line); ad = 2T and a g = OAgo (dashed 
line); ad = 4F and a g = 0.2go (dotted line); a d = 7r and 
a g — 0.2go (dashed dot line). 



could be used, for example, to transfer optically encoded 
digital information into a spin current. In a future publi- 
cation we plan to explore quantum noise induced switch- 
ing between stationary spin currents states for the single 
quantum dot case using a stochastic master equation^. 



VI. APPENDIX: LINEAR STABILITY 
ANALYSIS 



V. CONCLUSION 

In conclusion, we have studied the generation of spin 
currents by electrical transport through quantum dots 
coupled to a single mode of an optical microcavity. We 
have shown that when the cavity field is coherently 
driven, absorptive optical bistability occurs in the cavity 
field. Since the steady state spin current is proportional 
to both the cavity field amplitude and number of photons 
in the cavity, OB necessarily leads to bistability in the 
spin current. Moreover, we have shown that this bistabil- 
ity persists in the presence of inhomogeneous broadening 
of the dots' Zeeman splitting and vacuum Rabi frequen- 
cies. 

The cavity field and spin current could be made to 



In order to analyze the stability of the steady state 
solutions, we consider small fluctuations <5Pj, 5P\, Sa 
and 5s about steady the state solutions of Eqs. [T8ll2"Tl 

P 4 (t) = P^+SP^t) (32) 

P t (t) = P T (0) +<5P T (<) (33) 

a(t) = a (0) +Sa(t) (34) 

s(t) = s^+Ss(t). (35) 

where p[°\ P^° , a' ' and correspond to one of the 
three steady solutions solutions given by Eq. [5U By in- 
serting Eqs. I3"2"ll3"5l into Eqs. [TBll2"T1 and discarding terms 
that are quadratic in the fluctuations, we obtain the fol- 
lowing linear equations for the fluctuations, 
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The associated characteristic equation for the eigenvalues 
of J is a sixth order polynomial 

6 

5>,A Z =0 (37) 
i=i 

A particular steady state solution, p|°\ a' ' , s' ^}, 

is stable if Re[X] < for all six of the eigenvalues since 
any small noise induced fluctuation will then decay away 



exponentially. Steady state solutions with at least one 
eigenvalue satisfying Re[\] > are unstable since small 
fluctuations will grow exponentially with time. Numer- 
ical solutions for the eigenvalues for each of the three 
steady states in Figs. 4-6 indicate that two of the roots 
are stable and one is unstable. The stable roots corre- 
spond to the upper and lower branches in Figs. 4-6 that 
have positive slope while the middle branch with negative 
slope is unstable. 
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